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ON BACKSCATTER FROM SPATIALLY VARYING SURFACES 

Normal C. Grody 

ABSTRACT 

This paper deals with the theory of radar backscatter from 
rough (spatially varying) surfaces. An integral equation developed 
by Kerr to describe the radar backscatter from a perfectly conduct- 
ing surface is applied to non-smooth surfaces. The results are 
compared with those obtained by Beckmann as well as those obtained 
by Wright. Differences between the three solutions are discussed. 
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ON BACKSCATTER FROM SPATIALLY VARYING SURFACES 


INTRODUCTION 

The current interest in remote sensing of the environment has focused at- 
tention on the interaction of radiation with rough land and ocean surfaces. In 

recent years a number of theories have been developed with different approxima- 

1-9 

tions to determine the scattering from rough surfaces . Of the various ap- 

3 2 3 4 

proaches, the work of Beckmann etal ’ and Wright ’ have been most success- 
ful in explaining the radar backscatter measurements from sea and land 
surfaces^’ respectively. 

2 

Beckmann applied the Kirchhoff method using the tangent plane approxima- 

3 

tion to obtain the radar cross section, while Wright employed a zeroth order per- 
turbation technique. An alternate approach is presented here for a perfectly 
conducting surface which agrees basically with Beckmann's result, however a 
more accurate radar cross section is obtained by adjusting the end term contribu- 
tions to be exactly zero. Also, the small perturbation solution exhibits different 
angular dependence than Wright's result, although both methods agree under 
specific conditions. 


THEORY 


The relationship between transmitted power P t , and received power P r , of 

a radar system is given by the equation 

P r G 2 X 2 

— = o 

P t (47 r) 3 R 4 


( 1 ) 


1 



where G is the antenna gain, X the operating wavelength, and R is the distance 
between receiver to target. The quantity a is called the radar cross section and 
is given by the equation 



where k - 2i r /X . 

12 

These two equations were developed by Kerr for a uniform plane wave incident 
on a perfect conducting surface placed in a vacuum environment. Equation (2) is 
written for either horizontal or vertical polarizations, where effects due to other 
polarizations enter in the form of a multiplicative factor as given in Reference 
12. The z axes (refer to Fig. 1) corresponds to the direction of transmitted and 
received power and the differential area dA is that of the scattering surface pro- 
jected on a plane transverse to the z axes. 

The radar cross section was derived using the "tangent plane" approxima- 
tion, which physically assumes small wavelengths compared with the radius of 
curvature of the scattering surface. Hence, scattering from sharp discontinui- 
ties are not adequately treated using this formulation. Also, multiple scattering 
and shadowing effects are neglected in this approximation. These latter effects 
are especially important near grazing incidence so that the theory developed 
here does not apply for angles approaching grazing. 

Equation (2) is applied to the case of a rough conducting surface resulting 
from perturbations of a flat surface. The geometry used is shown in Figure 2, 
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Figure 2. Backscatter from a Perfectly Conducting Rough Surface 
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where the x, y, z coordinate system is defined by the direction of transmitted- 
received radiation (z axes) and the transverse plane xy on which the scattering 
surface is projected. The second set of coordinates x', y', z’ is defined such 
that the x'y' plane contains the unperturbed conducting surface whose normal 
direction is z'. 


Let 7 (x\ y') be the equation of the surface height at the location x', y'. 


The differential surface area dA is computed as follows: 


, , V' [z '-7 (x', y)] 

dA = a 7 • n dS , n = 

" ' IV' [z'-7(x',y)] 


so that. 


n = 


, 87 , fry 

8x' ~ y 8y' 


•§z - —x < il 


T + 


87 

8x 




also. 



where for small slopes, i.e. , 



dA 


( 97 

l cos0 , sin0 


9 y' 


) 


dx' dy' 


( 3 ) 
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It is observed that equation (3) is exact for one dimensional surface variations. 
Substituting (3) into (2), together with z = 7 (x\ y ') cos 9 + y' sin0, we obtain 



Equation (4) is the radar cross section for a rough surface whose height var- 
iation is given by the function 7 (x', y'). The function 7 is a measure of the per- 
turbation from a flat surface and may be considered a deterministic quantity (eg; 
sinusoidal reflecting grating) or a random variable (eg; sea surface). 


For large reflecting surfaces in comparison to the antenna coverage, the 
integrand in (4) must contain the gain function of the antenna. For simplicity, 
we will consider uniform coverage for y' = ±L/2, x' = ±L/2 and zero coverage 
outside this range. Hence, the integration limits are reduced as follows: 



where this result is simplified by integrating by parts the term 



d7 -\2k (y'sin 0 + 7 cos 0 ) 

ay' 


dx'dy' 


( 6 ) 
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..which can be written as. 


tan0 



(1 _ e -i 2 K 7 <:os 0 ) e -i2xy'sin0 dx'dy' 


(7) 


+ 



sin (fey cos0) 
fey cos0 


) 


e -i2« (y'sin0 + 'Ay cos 0) 


y'= L/2 
dx' 


J y'= -L/2 


where the above result was obtained by assuming y can be expressed as 
7 = f(x')g(y'). If we also assume g(± L/2) = 0, then the end term contributions 
are identically zero (for all angles) and equation (5) simplifies to 


a 


4n 

X 7 



— i2K (y'sin0 + y cos0) dx'^y' 


e -i2/cy' sin0 dx - dy - 


( 8 ) 


Beckmann et al^ applied the Kirchhoff method using the tangent plane approxi- 
mation to obtain the scattering from rough surfaces. By refering to Beckmann's 
solution2 of the radar cross section, it is found that the second integral term in 
(8) is missing. The reason for the missing term in Beckmann's result arises 
from the manner in which he evaluated the integral equation (6). His result. 
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obtained by integration by parts, is 

• L/2 


- tan0 ff e' i2K < y ' sin0 + 7 cos0) dx'dy' 

JJ- L/2 


(9) 


L/2 g -i2K (y'sinS + 7 cosy) 


r L/2 
J-L/2 


-i 2 k cos0 


y'= L/2 
dx' 

y'= -L/2 


where the end point contributions were simply neglected in comparison with the 
first term for kL >1. In order to obtain a more accurate result we utilized 
equation (7) which can be deduced from (9) by adding and subtracting the term 

r L l 2 . , . 

tan0 I e" l2Ky sin0 dx'dy' 

J- L/2 

As mentioned previously, equation (7) unlike (9) has zero end- term contributions 
when 7 = 0 at y' = ±L/2. 


The importance of the second integral term in (8) is easily seen by consider- 
ing a flat surface (7 = 0) so that we obtain the familiar result 


a 


4n 

X 2 


COS0 



2 


e ' 2Ky dx'dy' 



( 10 ) 


where A = L 2 is the area illuminated by the antenna and u = /cLsin0 . Had the 
second term been eliminated in (8), i.e. , Beckmann's result, the cos 2 0 factor 
would be replaced by sec 2 0 , resulting in an erroneous solution. 
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To determine the effects due to perturbations of a flat surface, the exponen- 
tial term in (8) containing 7 is expanded as follows: 


e -i2«7 cose 



(-i2<ccos9) n n 
n! 7 


(10 


which converges for the condition, 2kj cosd < 1. Substituting (11) into (8) we 
find 



00 

rr L /2 

2 a n 

| j 7 n e~ i2K y sin ^ dx'dy' 

n=0 J 

J -L /2 


( 12 ) 


where „ , (-i2xcos0) n 

a n = cost? and a n -A n = ; — . 

0 n ^° n! cos0 

This final form of a appears as a series of Fourier transforms of 7°. The re- 
strictions, as mentioned previously, are summarized below, 

(a) Tangent Plane Approximation 

(b) 7 = f (x') g(y') where g(±L/2) = 0 (13) 

(c) 2/<7 cosd < 1 

Equation (12) is useful for both small and large perturbations as long as the con- 
ditions (13) are met. The last part of this paper will consider small perturba- 
tions and compare the resulting equation with that obtained by Wright. 


SMALL PERTURBATIONS 

For small perturbations, i.e. , 2kj cos 0 « 1, keeping only the first two 
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terms of the series expansion (12) we obtain 


Q 

II 

k 


rr L/2 

rr L/2 

cos 9 

1 e ~ i2l< y sin0 dx'dy' - i2k 

7e -i2*y sin0 dx ’ dy ' 

X 2 

J 

J-L/2 J 

J -L/2 


(14) 

The first term is the contribution due to the unperturbed (flat) surface (Eq. 10) 
and is negligable for the condition 

|±0 1 > sin" 1 X/2L (15) 

where the right half of the inequality corresponds to the first zero crossing of 
equation (10). Hence, for angles sufficiently far from the zenith direction, the 
second term of (14) predominates so that we can write approximately 


4k 4 ff L/2 

"‘-■JJ.U2 


7e -i2Ky'sin0 dx ' dy ' 


L/2 

3 

This equation is similar to Wright's result when applied to a perfectly conduct- 
ing surface, except for the cos 4 0 multiplicative factor in Wright's solution. 

This difference can only be attributed to the different approximation made in 
Wright's analysis, i.e. , the use of the zeroth order perturbation solution for the 
scattered field rather than the direct use of the tangent plane- approximation. It 
is interesting to note that the cos 4 0 term can be obtained starting with equation 
(5) by retaining only the first order term cos0, disregarding the derivative term 
07 / 9 y'sin 0 and repeating the above analysis. 
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CONCLUSIONS 


The simple theory developed is useful for analyzing both small and large 
perturbations of a flat surface, under the restrictions enumerated above. It is 
found that this result is more accurate in describing the backscatter from rough 
surfaces than Beckmann’s solution. Upon comparing the small perturbation re- 
sult with Wright's formulation it is seen that they generally display different an- 
gular characteristics, although both agree under certain conditions. 

REFERENCES 

1. P. Beckmann and A. Spizzichino, "The Scattering of Electromagnetic Waves 
from Rough Surfaces", London/New York: Pergamon/Macmillian, 1963. 

2. P. Beckmann, "Scattering by Composite Rough Surfaces", Proc. IEEE, 

Aug. 1965, pp 1012-1015. 

3. J. W. Wright, "Backscatter from Capillary Waves with Application to Sea 
Clutter", IEEE Trans, on Antennas and Propagation, vol. AP-14, Nov. 1966, 
pp 749-754. 

4. J. W. Wright, "A New Model for Sea Clutter", IEEE Trans, on Antennas 
and Propagation, vol. AP-16, March 1968, pp 217-223. 

5. S. O. Rice, "Reflections of Electromagnetic Waves from Slightly Rough 
Surfaces", Commun. Pure Appl. Math. 4, no. 2/3, 1951, pp 351-378. 


10 



6. G. R. Valenzuela, "Scattering of Electromagnetic Waves from a Tilted 
Slightly Rough Surface", Radio Sci. , vol. 3, no. 11, Nov. 1968, pp 1057- 
1066. 

7. D. J. Fang, "Scattering from a Perfectly Conducting Sinusoidal Surface", 
IEEE Trans, on Antennas and Propagation, vol. AP-20, May 1972, pp 388- 
390. 

8. P. J. Lynch and R. J. Wagner, "Rough -Surface Scattering: Shadowing, 
Multiple Scatter, and Energy Conversion", J. Math. Phys. , vol. 11, no. 

10, Oct. 1970, pp 3032-3042. 

9. M. I. Sancer, "Shadow-Corrected Electromagnetic Scattering from a Ran- 
domly Rough Surface", IEEE Trans, on Antennas and Propagation, vol. 
AP-17, Sept. 1969, pp 577-585. 

10. N. W. Gunnard and J. C. Daley, "An Experimental Study of a Sea Clutter 
Model", Proc. IEEE, vol. 58, no. 4, April 1970, pp 543-550. 

11. P. Beckmann, "Radar Backscatter from the Surface of the Moon", J. 
Geophys. Res. , vol. 70, May 1965, pp 2345-2350. 

12. D. E. Kerr, "Propagation of Short Radio Waves", New York; McGraw-Hill, 
1951, p 697. 


NASA— GSFC 


11 



